FSU-SCRI-97-117 
September 1997 



1 



The Schrodinger Functional and Non-Perturbative Improvement^ 

R.G. Edwards a , U.M. Heller a and T.R. Klassen a 

a SCRI, Florida State University, Tallahassee, FL 32306-4130, USA 

After describing the Schrodinger functional for standard and improved gluon and quark actions we present 
results for the non-perturbative clover coefficients of the SW quark action coupled to the Wilson plaquette 
action for (5 > 5.7, as well as the Liischer-Weisz one-loop tadpole improved gauge action, both in the quenched 
approximation. 



1. Introduction 

The high cost of lattice QCD simulations has 
revitalized interest in the (on-shcll) improvement 
program. Within the Symanzik [Q approach, 
which we will follow, the use of one-loop (or even 
classical) and tadpole Q improved gauge actions 
has lead to much smaller scaling violations on 
coarse lattices than for the standard Wilson pla- 
quette action. Numerous studies of the static 
potential, thermodynamics, heavy quarks in ei- 
ther relativistic or non-relativistic frameworks, 
and glueballs (the latter on anisotropic lattices) 
have demonstrated this. References can be found 
in the LATTICE proceedings of the last few years. 

The improvement of quark actions is much 
harder. For Wilson-type quark actions, which 
we will consider here, this is ultimately due to 
the doubler problem. At least at the quantum 
level one incurs 0(a) violations of chiral symme- 
try, which have turned out to be quite large. 

A great step forward was recently taken 
by the ALPHA collaboration ||, which used 
the Schrodinger functional and the demand 
that the PCAC relation hold at small quark 
masses, to eliminate all on-shell 0(a) errors for 
Shcikholeslami-Wohlert (SW) f| quarks coupled 
to the Wilson gauge action. Various renormaliza- 
tion constants of axial and vector currents were 
also calculated non-perturbatively. 

The success of improved gauge actions on 
coarse lattices has motivated us to consider the 
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non-perturbative 0(a) improvement of quark ac- 
tions coupled to improved gauge actions. In the 
process we have also reconsidered the case of the 
SW action coupled to the Wilson gauge action 
and extended the determination of the 0(a) co- 
efficient to coarser lattices than in (H . 

Although we are also in the process of deter- 
mining the improvement coefficients of various 
currents, we will here concentrate on the 0(a) 
improvement coefficient of the action. Details of 
the general theoretical setup and our motivation 
can be found in [El; our results will be described 
in detail in future publications |@]. 

2. 0(a) and 0(a 2 ) Improvement 

For Wilson-type quark actions we have to in- 
troduce second order derivative and clover terms 
to eliminate doublers without introducing classi- 
cal 0(a) errors. On the quantum level, on an 
isotropic lattice, these two terms are still the 
only ones that exist at 0(a). We write them as 
ar(J2^ + ^iva-F). One of the coefficients r, u> 
can be adjusted at will by a field transformation. 
It is convenient to fix the Wilson parameter r; to 
eliminate all 0(a) violations of chiral symmetry 
we then have to tunc the clover coefficient u as a 
function of the gauge coupling. 

Note that the 0(a) terms in the action break 
chiral but not rotational symmetry (at this or- 
der), whereas the leading 0(a 2 ) errors, that al- 
ready exist at the classical level, show the op- 
posite behavior; they break rotational but not 
chiral symmetry. For this reason the 0(a) and 
leading 0(a 2 ) terms can essentially be tuned in- 
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dependency (cf. ) . By the same token, one can 
argue that one indeed should tune the 0(a) and 
(leading) 0(a 2 ) terms to eliminate the violations 
of both chiral and rotational symmetry. 

Eliminating the leading 0(a 2 ) errors in a quark 
action leads to the D234 actions || . As for gauge 
actions, using classical and tadpole improvement 
at 0(a 2 ) seems to almost completely eliminate 
the violation of rotational symmetry ||. 

So far we have discussed isotropic lattices. 
Anisotropic lattices, with a smaller temporal 
than spatial lattice spacing, are of great inter- 
est for studies of heavy particles (glueballs, heavy 
quarks, hybrids) and thermodynamics. Improve- 
ment is more complicated for actions on such lat- 
tices. After considering the most general field 
redefinitions up to 0(a), one sees [|J that two 
more parameters have to be tuned for on-shell 
improvement of a quark action up to 0(a). One 
already appears at 0(a°), namely, a "bare ve- 
locity of light" that has to be tuned to restore 
space-time exchange symmetry (by, say, demand- 
ing that the pion have a relativistic dispersion re- 
lation for small masses and momenta) . The other 
is at 0(a); the two terms that have to be tuned 
at this order can be chosen to be the temporal 
and spatial parts of the clover term. 

Although the general methods sketched here 
should eventually be useful also for the 
anisotropic case, we will in the following restrict 
ourselves to isotropic lattices. 

3. Chiral Symmetry Restoration 

Consider QCD with (at least) two flavors of 
mass-degenerate quarks. The idea || for deter- 
mining the clover coefficient is that chiral sym- 
metry will hold only if its Ward identity is sat- 
isfied as a local operator equation. In Euclidean 
space this means that the PCAC relation between 
the iso- vector axial current and the pseudo-scalar 
density, 

(d,Al(x)0)=2m(P b (x)0) (1) 

should hold for all operators O, global boundary 
conditions, x (as long as x is not in the support 
of O) etc. More precisely, it should hold with 
the same mass m up to 0(a 2 ) errors (which are 



quantum errors for (classically and tapole) a 2 im- 
proved actions). This will only be the case for the 
correct value of the clover coefficient lo. 

Several issues have to be addressed before this 
idea can be implemented in practice. First of all, 
even though here we can ignore the multiplicative 
renormalization of A b and P b , there is an additive 
correction to A b ^ at 0(a), 

P b (x) cx i>(x) l5 ±T b iP(x) , 

A\{x) cx ^(x) 1 ^T b ^(x)+ac A d tl P b (x)(2) 

The determination of lo is therefore tied in with 
that of ca- We will see later how to handle this. 

Note that lo and ca have an 0(a) ambiguity 
(at least at the quantum level); different improve- 
ment conditions will give somewhat different val- 
ues for lo and ca ■ Instead of assigning an error to 
lo and ca one should choose a specific, "reason- 
able" improvement condition — the associated 
0(a 2 ) errors in observables are guaranteed to ex- 
trapolate away in the continuum limit. 

For various conceptual reasons it is preferrable 
to impose the PCAC relation at zero quark mass. 
Due to zero modes this is not possible with peri- 
odic boundary conditions (BCs); the quark prop- 
agator would diverge. Another reason to abandon 
periodic BCs is that to be sensitive to the value 
of lo it would be highly advantageous to have a 
background field present; it couples directly to the 
clover term. The Schrodinger functional provides 
a natural setting to implement these goals. 

4. The Schrodinger Functional 

The phrase "Schrodinger Functional" (SF) 
refers to quantum field theory with Dirichlet, 
i.e. fixed, BCs ||. In the following we will al- 
ways use periodic BCs in space (extent L) and 
fixed BCs in time (extent T). For finite T the 
Dirac operator has a gap of order 1/T even for 
vanishing quark mass, at least at weak coupling. 
Furthermore, by choosing different BCs at "oppo- 
site ends of the universe" one induces a chromo- 
electric classical background field. 

In implementing the SF on the lattice the main 
point is to understand exactly how to impose 
fixed BCs on the gauge and quark fields. In par- 
ticular, we must be able to do so for improved 
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Figure 1. Temporal weight factors of various pla- 
quettes in the Wilson gauge action on a lattice 
with T = 5a and L = 5a. 



actions. For details we refer to S; here we just 
mention some salient features: 

1. The main difference between the Wilson 
and 0(a 2 ) improved gauge actions is that for the 
latter the "boundary" consists of a double layer 
of time slices. To avoid boundary errors larger 
than those of the bulk action, one must, already 
at the classical level, assign loops at the boundary 
special "temporal weight factors" that depend on 
the temporal extent of the loop (cf. figs. 1 and 2). 
The classical values of the weight factors are easy 
to understand from elementary calculus formulas 
(e.g. the trapezoidal rule explains the factors of \ 
in the Wilson case of fig. 1). When using the SF 
as a tool to tune coefficients in a local action (or 
current), it is fortunately not necessary to know 
the exact quantum values of the boundary coef- 
ficients: the local Ward identities have to hold 
independent of global effects at the boundary. 

2. If the boundary values of the gauge field at 
the top and bottom of the universe commute, the 
following is a solution of the lattice field equations 
for any gauge action: Uq(x) = 1, and 



' "/. (•'■) = <"M> ( j; [xoC' k + (T- xo)C k 



(3) 



(The boundary values of the gauge field can be 
read off from the above by evaluating it at xq = 
and T for Wilson, respectively, xq = 0, a and 
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Figure 2. As in fig. 1 for an improved gauge ac- 
tion. The thick lines indicate the inner boundary 
layers. To avoid classical boundary errors one 
must choose b = 1/24. 



T,T + a in the improved case.) 

The question is if the above background field U 
is the unique (up to gauge equivalence) absolute 
minimum of the classical action for given bound- 
ary values. Uniqueness is important, e.g. for per- 
turbative calculations. A theorem establishing 
uniqueness holds in the Wilson case ||, if the 
Cfc, C' k parameterizing the boundary values sat- 
isfy certain conditions. In the improved case it 
has been checked using simulated annealing that 
uniqueness holds under the same conditions 

3. To impose consistent fixed boundary condi- 
tions for the fermion fields, it is sufficient to con- 
sider the projector structure of the field equations 
(more precisely, it is only the projector structure 
in the time direction that matters). For an action 
with the same projector structure as the standard 
Wilson quark action, one has to specify P + ip{x) : 
ip(x)P~ at the (inner) lower boundary in figs. 1 
and 2, and P~ip(x), ip(x)P + at the (inner) up- 
per boundary. Here P^ = i(l ± 70). For an 
improved quark action with the appropriate pro- 
jector structure || one has to specify the same 
components on both the inner and outer bound- 
ary layers in fig. 2. 

4. One of the very useful ideas in applications 
of the SF is that of quark boundary fields B . They 
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are defined as functional derivatives, within the 
path integral, with respect to the boundary val- 
ues specified in the previous paragraph (which are 
then set to zero). For the improved case one can 
actually define two sets of quark boundary fields. 
It turns out that if one defines them with respect 
to the outer boundary values, then most formulas 
relevant for our application of the SF are identical 
for improved and standard actions. The bound- 
ary fields corresponding to the above boundary 
values will be denoted as C( x )> C( x ); C'( x ); C'( x ); 
the first pair being the lower, the second the up- 
per boundary fields. 

5. Details of Non-Perturbative Tuning 

With O b = a 6 J2 yz CO) 75 5A(z), in terms of 
the lower boundary fields, and 
1 



the PCAC relation becomes 

m(xo) = r(x ) + acA s(x ) = m + 0(a 2 ) . (5) 

where rh differs from m by irrelevant multiplica- 
tive renormalization factors, and 



r(sco) 



s(x ) 



(G) 



2/ P (xo) v 2/ P (xo) 
Here Vo and Ao are standard first and second 
order lattice derivatives, respectively (in the im- 
proved case one actually has to use an improved 
first order derivative to be consistent). 

Similarly, f' Xl m / ,r',s' are defined in terms of 
the upper boundary fields C ■ 

From m(yo) = m'(zo) one obtains an estima- 
tor of ca '■ 



- / \ 1 r(y )-r'(z ) 

CA(yo,z ) = — 



(7) 



a s{yo) - s'(zo) ' 
In terms of a suitable ca we now have two differ- 
ent estimates of the current quark mass: 

M(x Q ) = r(x ) + a c A s(x ) , 
M'(x ) = r'Oo) + ac A s'(x ). (8) 

Their equality in the presence of a suitable back- 
ground field will be our improvement condi- 
tion for to. More precisely, we demand 

AI(i ) = M(x ) - M'{x a ) = AM (0) (x ) (9) 



for some well-chosen xo; here the superscript (0) 
denotes the higher order (and small) tree-level 
value of the quantity in question. In practice one 
measures the required correlators in a simulation 
for several trial values of uj, and interpolates to 
find the zero crossing of AM — AM' '. This de- 
termines the non-perturbative value of u>. 

6. Results 

The results we describe in the following all refer 
to the SW action on cither Wilson or one-loop 
tadpole improved glue |l(],[ll| (which we will refer 
to as "LW glue"). We will always work in the 
quenched approximation on isotropic lattices. 

Below we use some CA(yo,yo) as estimator of 
ca in eq. (§). We then denote AM(x ) by 
AM(x Q , y ). We used lattices T ■ L 3 = 15 • 8 3 or 
12 -6 3 for Wilson glue, and 14 -8 3 for LW glue. Af- 
ter some study we decided to use AM(12, 4) and 
AM (9, 3), respectively, in the improvement con- 
dition for lu in the Wilson case. In the improved 
case we chose AM(11, 5). Typically we generated 
1000—2500 configurations for each gauge coupling- 
considered. 

Even though the SF alleviates problems due to 
zero modes, it turns out that on coarse lattices 
fluctuations still lead to accidental zero modes 
at vanishing quark mass ( "exceptional configura- 
tions"). Fortunately, it turns out that the mass 
dependence of the non-perturbative u is so weak 
that one can safely determine it at larger mass 
values. In this manner we have extended the non- 
perturbative clover coefficient obtained by the 
ALPHA collaboration for > 6.0 to j3 > 5.7. 
In fig. 3 we illustrate the weak mass (and vol- 
ume) dependence of u for (3 — 5.7 Wilson glue. 
We have also checked that the mass dependence 
is weak for (3 = 5.85 and 6.0. The same can be seen 
for LW glue in fig. 4, which furthermore demon- 
strates the linearity of AM as a function of u> (all 
results shown in fig. 4 were calculated using the 
same gauge configurations). 

For future use of improved quark actions it is 
advisable to present the non-perturbative clover 
coefficient as a definite function of the bare gauge 
coupling. Combining our Wilson results for [3 = 
5.7, 5.85 and 6.0 with those of the ALPHA col- 
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Figure 3. The non-perturbative clover coefficient 
co for (3 = 5.7 Wilson glue for different quark 
masses and volumes, and its extrapolation to 
M=0. 
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laboration and one-loop perturbation theory we 
obtain the parameterization (g 2 = 6//3) 



,. 1 - 0.60503 2 - 0.1684 5 4 



u(9 2 ) = 



1 - 0.8709 5 2 



g 2 < 1.06(10) 



We were able to accommodate our (3 — 5.7, 5.85 
clover coefficients and the value for j3 — 6.2 
from in our curve only in a slightly unsatis- 
factory manner (extending the Pade in either the 
numerator or denominator does not help). This 
issue is under investigation. In the interim the 
curve (|Io| ) should be regarded as preliminary. 

For the case of LW glue our current data are 
parameterized well by 



^lw(5 2 



1 - 0.3590 g 2 
1 - 0.4784 g 2 



g 2 < 1.55. 



(11) 



(The one-loop coefficient is presently not known 
analytically in this case.) The relation between 
the coefficient of the plaquette term in the LW 
action, /3lw, and the bare coupling g 2 is, to one 
loop JH}, g 2 = 10/(/3 LW - 1.422). For larger cou- 
plings than those in ([n]) it seems that the non- 
perturbative LOhw determined with the SF rises 
dramatically. This is currently under investiga- 
tion. 

It is interesting to compare the non- 
perturbative clover coefficients obtained for Wil- 
son and LW glue at the same physical scale. The 
string tension has been measured for both of these 
actions, so we will use it to set the scale. We 
would like a curve parameterizing the string ten- 
sion a as a function of the coupling. It is known 
that the two (or three) loop running of the cou- 
pling does not properly describe it's lattice spac- 
ing dependence (nor that of other observables) for 
the couplings we are interested in. However, as 
pointed out in |Q, this is not to be expected, 
since a has discretization errors, g 2n a 2 , a 4 , . . . 
(n = 0/2 for Wilson/LW glue) that should be 
taken into account. We therefore try to parame- 
terize (a^fa){g) as 



Figure 4. Results from the non-perturbative 
tuning of to at different quark masses for LW glue 
at /3lw=8.4 (where a « 0.1 fm). 



X fi9 2 ) ( 1 + c 2 9 2n a(g) 2 + c 4 a(.g) 4 ) , (12) 

in terms of the three fit parameters y/a/A, C2, C4. 
Here a(g) = f{g 2 )/f({), in terms of the universal 
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Figure 5. Non-perturbative clover coefficients 
versus lattice spacing; scale set by y/a — 460 MeV. 
The dashed line denotes the curve from || . Small 
symbols indicate the tree-level (plaquette) tad- 
pole estimates for w. 



two-loop function 

a{g)K = f{g 2 ) ee (b g 2 )'^ exp(- 



-) (13) 



This works very well (for details see 0), and 
the clover coefficients arc presented as a function 
of the lattice spacing in fig. 5. It is interesting to 
observe that for a > 0.05 fm both the tadpole and 
the non-perturbative a/s are pretty much linear 
in a, at least up to about 0.15 fm. Furthermore, 
the differences between the non-perturbative and 
the tadpole values are essentially linear down to 
the smallest couplings considered for both cases 
(of course, for sufficiently small couplings the dif- 
ferences should be of order g 2 ). 

7. Conclusions and Outlook 

We have shown that the Schrodinger functional 
and the non-perturbative elimination of O(a) er- 
rors in a Wilson-type quark action can be suc- 
cessfully extended to improved (gauge) actions. 
By establishing that the non-perturbative clover 



coefficient has a very weak mass dependence we 
were able to determine it for lattice spacings sig- 
nificantly above 0.1 fm, for both Wilson and im- 
proved gauge actions. 

We are currently investigating different defini- 
tions of the axial current improvement coefficient 
ca ■ This is necessary if one wants to determine it 
on coarse lattices, where some definitions lead to a 
ca of rapidly increasing magnitude (at least when 
using Wilson glue) . Once the ca determination is 
completed we plan to determine the other current 
normalization and improvement coefficients ||. 

There are many other situations in which the 
non-perturbative elimination of the O(a) viola- 
tions of chiral symmetry is important, the most 
obvious examples being full QCD and D234 
quarks on anisotropic lattices (for the study of 
heavy quarks). We hope that the ultimate out- 
come of this and future studies will be the abil- 
ity to perform accurate continuum extrapolations 
from much coarser, and therefore cheaper, lattices 
than hitherto possible. 
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